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relative  contributions  of  surface-catalyzed  and  gas-phase  reactions  to  the  total 
heat  release  race  in  a  catalytic  combustor  over  a  range  of  operating  variables 
that  include  gas  flow  velocity,  catalyst  temperature,  surface  reactivity,  and 
fuel-air  ratio.  Two  catalytic  combustor  systems  were  employed.  One  involved 
stagnation  point  boundary  layer  flow  with  a  catalytically  active  flat  surface, 
the  other  duct  flow  in  a  multichannel  catalytic  monolith.  The  experimental 
studies  were  supplemented  by  theoretical  calculations  of  heat  release  rates  and 
temperature  profiles  based  on  models  appropriate  to  the  geometry  and  fluid  dynamips 
of  each  catalytic  combustion  system.  The  result  of  our  studies  with  dilute 
propane-air  gas  mixture  and  platinum  coated  surfaces  have  demonstrated  that  in 
both  combustor  configurations  the  gas  phase  combustion  can  be  Initiated  by  the 
wall-catalyzed  react ion. stagnation-point  flow  the  heat  release  due  to  the 
homogeneous  reaction  is  significant  at  low  gas  velocities  and  high  surface 
temperatures.  In  this  rggusn  of  predominant  gas-phase  reaction  the  temperature 
profile  exhibits  a  maximum  wi\h  respect  to  distance  from  the  catalyst  surface. 

In  the  miltichannel  monolith  combustor,  the  axial  temperature  profile  inside  a 
centrally-located  duct  is  strong^  affected  by  heat  conduction.  Results  of  the 
computational  analysis  provide  valuable  information  on  the  effect  of  various 
parameters  on  the  temperature  and  reaction  distribution  ir.  a  catalytic  monolith 
Including  surface  reaction  rate,  gas-phase  reaction  rate  transport  parameters, 
axial  conductivity,  and  flow  velocity .  \The  temperature-profile  characteristic 
observed  experimentally  correspond  satisfactorily  to  those  predicted  by  the 
model  for  the  case  of  exothermic  heat  release  by  surface  initiated  gas-phase 
reaction.  The  report  is  divided  into  thrde  sections:  (1)  "Catalytic 
Combustion  in  a  Stagnation  Point  Boundary  Layhr"  (2)  "Combustion  in  a  Monolithic] 
Catalytic  Combustor  with  Gas  Phase  Reaction"  (-3)  "Numerical  Analysis  of  Catalytip 
Monolith  Combustion" 
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SUMMARY 

Engineering  studies  with  catalytic  combustors  in  the  main  burner 
and  afterburner  of  aircraft  propulsion  systems  have  been  concerned  with 
the  attainment  of  high  combustion  efficiency  and  flame  stabilization  of 
fuel-air  mixtures  outside  the  conventional  flammability  limits.  The 
development  of  such  devices  that  will  perform  at  optimum  efficiency 
over  a  broad  range  of  operating  conditions  requires  a  detailed  under¬ 
standing  of  the  flow  of  a  reactive  gas  through  a  duct  with  catalytic 
walls.  Both  heterogeneous  (wall)  and  homogeneous  (gas)  chemical  reaction 
kinetics  as  well  as  the  fluid  dynamics  of  laminar  and  turbulent  flow 
need  to  be  included  in  a  model  of  the  catalytic  combustion  system.  To 
apply  such  a  model  it  is  necessary  to  establish  the  chemical  kinetic 
parameters  for  the  surface  catalyzed  and  the  gas-phase  reactions  that 
individually  or  jointly  govern  the  heat  release  rates.  Of  special 
interest  to  our  study  is  the  evaluation  of  the  relative  contributions 
of  surface-catalyzed  and  gas-phase  reactions  to  the  total  heat  release 
rate  in  a  catalytic  combustor  over  a  range  of  operating  variables  that 
include  gas  flow  velocity,  catalyst  temperccure,  surface  reactivity, 
and  fuel-air  ratio. 

In  order  to  examine  the  respective  roles  of  the  various  operating 
parameters  in  the  combustion  process,  we  studied  two  catalytic  combustor 
systems.  One  involved  stagnation  point  boundary  layer  flow  with  a 
catalytically  active  flat  surface,  the  other  duct  flow  in  a  multichannel 
catalytic  monolith.  The  experimental  studies  were  supplemented  by 
theoretical  calculations  of  heat  release  rates  and  temperature  profiles 
based  on  models  appropriate  to  the  geometry  and  fluid  dynamics  of  each 
catalytic  combustion  system. 

The  result  of  our  studies  with  dilute  propane-air  gas  mixture  and 

platinum  coated  surfaces  have  demonstrated  that  in  both  combustor 

configurations  the  gas  phase  combustion  can  be  initiated  by  the  wall- 

catalyzed  reaction.  In  stagnation-point  flow  the  heat  release  due  to 
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the  homogeneous  reectlon  is  significent  at  low  gee  velocities  end  high 
surface  temperatures.  In  this  region  of  predominent  ges-phase  reectlon 
the  temperature  profile  exhibits  a  maximum  with  respect  to  distance  from 
the  catalyst  surfaca.  In  the  multichannel  monolith  combustor,  the  axial 
temperature  profile  inside  a  centrally-located  duct  is  strongly  affected 
by  heat  conduction.  Under  our  experimental  conditions  the  wall  reaction 
is  controlled  by  surfaca  kinetics  and  not  by  mass  transport  of  reactants 
to  the  wall.  The  location  of  the  maximum  temperature  and  its  magnitude 
is  a  function  of  the  mass  flow  rate  of  reactant  through  the  combustor. 
These  experimental  data  are  interpreted  with  the  aid  of  analytical 
modeling. 

Results  of  the  computational  analysis  provide  valuable  information 
on  the  effect  of  various  parameters  on  the  temperature  and  reactant 
distribution  in  a  catalytic  monolith, including  surface  reaction  rate, 
gas-p;..ase  reaction  rate,  transport  parameters,  axial  conductivity,  and 
flow  velocity.  The  temperature-profile  characteristics  observed  experi¬ 
mentally  correspond  satisfactorily  to  those  predicted  by  the  model  for  the 
case  of  exothermic  heat  release  by  surface  initiated  gas-phase  reaction. 

To  attain  high  conversion  efficiency  of  the  fuel  within  a  short  length 
of  crtalyst,  it  is  necessary  to  have  high  rates  of  surface-catalyzed 
and  gas-phase  reactions.  Increasing  the  reactivity  of  the  solid  surface 
does  not  necessarily  improve  performance  because  of  transition  to  a 
diffusion-limited  mode  of  operation.  Very  hi^a  mass  flux  of  the  fuel- 
air  mixture  is  detrimental  to  performance  because  it  shortens  the  resi¬ 
dence  time  and  reduces  the  contribution  of  the  homogeneous  reaction. 

The  following  sections  of  this  report  include  a  detailed  description 
of  the  experiments  and  the  results  of  our  study.  These  sections  are 
in  the  form  of  three  Technical  Reports  entitled: 

(1)  "Catalytic  Combustion  in  a  Stagnation  Point  Boundary 

Layer" 

(2)  "Combustion  in  a  Monolithic  Catalytic  Combustor  with 

Gas  Phase  Reaction" 

(3)  "Numerical  Analysis  of  Catalytic  Monolith  Combustion" 
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CATALYTIC  COMBUSTION  IN  A 
STAGNATION  POINT  BOUNDARY  LAYER 

by 

C.  M.  AbloV;  S.  Schechter,  end  H.  Wise 
SRI  International 
Menlo  Park,  California  94025 

ABSTRACT 

In  order  to  examine  the  relative  contribution  of  surface  and  gas 
phase  reactions  to  the  exothermic  conversion  of  a  fuel-air  gas  mixture 
flowing  over  a  catalytic  surface,  we  have  carried  out  a  theoretical 
and  experimental  study  of  catalytic  combustion  under  stagnation  point 
flow  conditions.  In  the  presence  of  exothermic  surface  reaction  the 
theoretical  model  provides  an  analytical  solution;  in  the  presence  of 
both  reaction  modes,  homogeneous  and  heterogeneous, computer  solutions 
were  obtained  for  the  total  heat  flux  at  the  surface  and  the  distribution 
of  temperature,  reactant,  and  product  concentrations  in  the  stagnation 
point  boundary  layer.  Having  available  kinetic  data  on  the  platinum- 
catalyzed  oxidation  of  propane  and  on  the  gas  phase  reaction  between 
propane  and  air,  we  selected  this  chemical  fuel  system  for  experimental 
and  theoretical  study.  Experimentally  the  flow  of  the  propane-air  mix¬ 
ture  (1  volZ  C^Hg)  was  directed  at  a  quartz  plate  whose  surface  was 
coated  with  thin  strips  of  vacuum-deposited  platinum,  that  served  both 
as  catalyst  and  resistance  thermometer.  During  an  experiment  the  temper¬ 
ature  of  each  strip  was  maintained  constant  (±  5  K)  by  adjusting  the 
electrical  heat  input  as  monitored  by  resistance  measurements.  By  this 
procedure  the  heat  released  by  exothermic  reaction  of  the  fuel-air  mixture 
could  be  determined  from  the  difference  in  electrical  power  required  to 
keep  each  Pt  stlp  at  its  original  temperature.  With  the  aid  of  the 
theoretical  analysis  we  were  able  to  compute  the  fractional  contributions 
of  catalytic  and  gas  phase  combustion  to  the  total  heat  flux  conducted 
to  the  catalytic  surface.  The  relative  contribution  of  each  combustion 
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nods  depended  on  cetelyet  activity,  volumetric  flow  rate,  end  fuel/eir 
ratio.  Under  our  experimental  conditions  with  transport  limited  surface 
reaction  an  Increase  in  reactant  flow  rate  enhancee  the  surface-catalysed 
contribution  to  the  total  heat  release  rate.  The  results  obtained  clearly 
indicate  the  conditions  under  which  either  combustion  mode  predominates. 


I  INTRODUCTION 


The  development  of  catalytic  reactors  for  efficient  combustion  uf 

various  fuel-air  mixtures  is  receiving  considerable  attention  not  only 

1  2 

in  pollution  control  but  also  in  gas  turbine  thruster  technology.  * 

To  a  large  measure  the  utility  of  catalytic  combustors  rests  upon  their 

ability  to  operate  efficiently  over  a  much  wider  range  of  fuel-air  ratios 

chan  are  imposed  by  the  flammability  limits  of  conventional  gas  buj-ners. 

More  recently  the  use  of  catalytic  devices  has  been  examined  for  NO^- 

free  combustion  in  gas  turbine  engines.  Another  interesting  application 

in  aircraft  propulsion  systems  is  flame  stabilization  in  after  burners 

3  4 

by  means  of  catalytic  devices.  ’ 

While  the  selection  of  catalytic  materials  for  efficient  combustion 
has  received  considerable  attention,^  an  analysis  of  the  relative  contri¬ 
bution  of  gas-phase  reactions  and  surface-catalyzed  processes  to  the 
performance  of  catalytic  combustors  over  a  range  of  operating  conditions 
is  made  difficult  by  the  complex  interplay  of  fluid  dynamics  and  reaction 
kinetics.  The  complexity  of  the  theoretical  analysis  is  lessened  some¬ 
what  for  the  case  of  catalytic  combustion  with  stagnation  point  flow 
conditions,  because  the  flow  field  applicable  to  such  a  system  is  well 
understood.®  The  reactive  gas  mixture  chosen  "or  our  studies  was  com¬ 
posed  of  propane-air  since  the  global  kinetics  for  the  platinum-catalyzed 

7  8 

surface  reaction  and  for  the  homogeneous  gas-phase  reaction  have  been 

measured  in  our  laboratory.  The  theoretical  model  for  reactive  stag¬ 
nation  point  flow  over  a  catalytic  surface,  the  experimental  measurements 
of  heat  release  rate,  and  the  interpretation  of  the  results  in  terms  of 
the  model  are  presented  in  the  following  sections. 
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II  THEORETICAL  MODEL 


In  eh*  theoretical  modeling  of  eh*  stagnation  point  flow  catalytic 
reactor  wa  consider  tvo-dimenaional,  steady  flow  of  a  compressible,  vis¬ 
cous  fluid  directed  perpendicularly  against  a  flat  solid  surface.  The 
fluid  is  composed  of  a  mixture  of  gaseous  fuel  and  ait  that  can  undergo 
exothermic  reaction  in  the  gas  phase  and  on  the  surface  of  the  solid 
by  catalytic  reaction. 

Cartesian  coordinates  (x,y)  are  Introduced  with  origin  at  the 
stagnation  point,  the  x-axls  along  the  surface  in  an  outflow  direction, 
and  the  y-axis  in  the  direction  of  the  oncoming  flow.  The  velocity 
components  are  u  and  v  in  the  x  and  y  directions,  respectively.  From 
the  symmetry  of  the  flow  field,  u  will  be  an  odd  function  of  x,  while 
the  other  dependent  variables,  v,  density  p,  pressure  p,  temperature  T, 
and  fuel  mass  fraction  Y,  will  be  even  functions  of  x.  Near  the  y-axis 

u  -  x  Ux(y)  +  ... 

v  -  VQ(y)  +  x2V2(y)  +  ...  (1) 

1 

p  -  pQ(y)  +  *  p2(y)  +  ••• 

where  the  subscript  is  the  same  as  the  power  of  x.  Subscripts  0  and  1 
have  been  omitted  below  for  brevity. 

The  equation  for  conservation  of  total  mass, 

d-M,  +  .  o 

dx  Sy 

becomes  to  the  lowest  order  in  x, 

pU  +  (pV)'  -  0  (2) 

where  the  prime  denotes  differentiation  with  respect  to  y.  The  Navier- 

9 

Stokes  equations  for  conservation  of  momentum  reduce  to 
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t 


(3) 


p(U2  +  VU')  -  OiU'y  +  (1/2)  p2  -  0 

pW'  -  (4/3)  GiV')'  +  p'  -  0  .  (4) 

Equation  (4)  shows  that  pressura  and  valocity  variations  go  togathar. 
For  tha  slow  flows  of  tha  praaant  experiments,  the  pressure  variation 
is  a  small  fraction  of  its  average  value  so  that  Equation  (4)  may  be 
replaced  by 


p  •  constant  (5) 

Variations  of  p2  with  y  are  similarly  negligible,  so  that  p^  in  Equation 
(3)  is  taken  to  be  a  constant. 

The  equation  for  conservation  of  fuel  reads ^ 

“(“t!  +  t?)  "a! (pl)!j)‘at(pDt?)‘  '*0  (6) 

where  D  is  the  diffusion  coefficient  and  Rg  is  the  rate  of  consumption 
of  fuel  by  gas-phase  reaction.  Near  the  solid  surface,  derivatives  with 
respect  to  x  are  negligibly  smaller  than  those  with  respect  to  y.  To 
the  lowest  order  in  x,  one  obtains 

pVY'  -  (p DY * Y  -  -RQ  .  (7) 

Fuel-lean  conditions  are  assumed  so  that  variations  in  the  concentrations 
of  oxygen  and  inert  species  need  not  be  considered. 

If  the  contributions  of  mechanical  work  and  viscous  dissipation  are 

9 

neglected,  energy  conservation  requires  for  a  steady  state  that  heat 
production  be  balanced  by  convection  and  diffusion. 

P  [« HP-  ♦ * *&]  -  si  f£)  -  h  W  •  • 

where  C  is  the  specific  heat  at  constant  pressure,  \  the  thermal 
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conductivity,  and  (AH)  tha  haat  of  exothermic  raaction.  Near  tha  stag¬ 
nation  point  ona  obtains 


pV(CT)'  -  (ll')'  -  (AH)  Rq  (9) 

Sines  tha  prassura  is  constant  by  Equation  (5)  tha  parfact  gas  law  givas 

pi  ■  constant  .  (10) 


At  y  *  »,the  variablas  have  tha  known  valuas 

T  -  T*  .  P  «  P.  ,  Y  «  ,  U  -  p  .  (11) 

Semlempirical  expressions'*'1  determine  3  for  given  nozzle  geometry  and  flow. 
Since  U/  •  *0,  Equation  (3)  gives 


p2  ■  -2p<y  •  <i2> 

On  the  surface  of  the  solid,  at  y  *  0,  tha  velocity  vanishes,  tha 
temperature  is  prescribed,  and  the  surface  reaction  is  fed  by  diffusion: 

U  -  v  -  0  ,  T  -  Ts  ,  D(pY)‘  -  Rs  ,  (13) 

where  Rg  is  the  rate  of  consumption  of  fuel  by  the  surface-catalyzed 
heterogeneous  reaction. 

Both  reaction  rates  are  taken  to  be  of  first  order  in  the  fuel 
concentration  and,  for  these  oxygen-rich  conditions,  of  zero  order  in 
oxygen ; 


Rq  "  PY  ,  Rg  *  kg  pY 
Nondimensionalization 

12 

Following  the  conventional  boundary  layer  development  we  introduce 
the  nondimens ional  independent  variable  T\  and  mass  flux  f,  a  function  of 
Tl,  by 
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f  -  -pv^p^p  JT 


(13) 


Than  by  Equation  (2) 

PU  "  P.»^  (16) 

whara  tha  subacrlpt  1}  danotaa  dlffarantiatlon.  Equation  (7)  becomes 

f  ^  +  a/**.)  <PD1tn>^  •  kG  pY/M  (17) 

If  tha  variation  of  C  la  neglected  ao  that 

CT'  -  (CT)'  ,  (18) 

Equation  (8)  may  ba  written 

fTq  +  (1  /V.J  (LepDX^  -  -kQ  (3H)pY/PoB8C  ,  (19) 

whara  La  is  tha  Lewis  number,  Le  ■  A. /p CD . 

Elimination  of  y,  U,  and  V  from  Equation  (3)  by  means  of  Equations 
(12),  (15)  and  (16)  gives 

^  -  ofoyp),,  -  (p/i*.)  c^yp^Ti  - 1  * 0  (20) 

Boundary  conditions  on  Equations  (17),  (19)  and  (20)  as  given  in 
Equations  (11)  and  (13)  become 

T  -  To  ,  p  -  p.  ,  Y-V,  ,  Ejj-1  at  1]«« 

1/2 

f  ■  fijj  ■  0  ,T«TS  ,  (pY)^  "  (Rs/D)  (^./P.6) 

at  T\  -  0  .  (21) 
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V  tha  fractional  variation  of  tamparatura  la  snail,  the  variations 
of  p  and  \k  In  Equation  (20)  can  ba  naglactad  to  obtain  an  aquation  for 
f  alona. 

*TO  +  tfWl  '  *i|  +  1  '  0  •  <22> 

This  aquation  with  tha  boundary  conditions  on  f  in  Equation  (21)  has 

6,12 

baan  solved  nunarlcally  by  Hlamanz  and  others.  Tha  tabulated  results 
show  that,  at  1]  ■  2.0,  and  u  have  reached  97Z  of  their  values  at 
T]  ■  The  boundary  layer  is  therefore  taken  to  extend  from  T)  ■  0  to 
H  •  2.0.  A  convenient  analytic  approximation  to  f  that  has  two  place 
accuracy  throughout  the  boundary  layer  is  given  by 

f  -  fal]2  +  bi]3)/(l  +  cTJ  +  bi]2) 

a  -  0.6047  ,  b  -  0.2692  ,  c  -  0.6338  .  (23) 

The  indefinite  integral  of  f  is  also  needed: 


2\ i  /ov 2 


j'f  dTl-Tlz/2  +  (a-c)T]/b  -  b  +  c(a-c)][ln(l  +  cTJ  +  W|  )]/2b 


+  £ab  -  c(3b  +  ac  -  c2)][arc  tan  (c  +  2bT| ) /s] /b2s  , 


...  2.1/2 
s  *  (4b  -  c  ) 


When  constant  average  values  are  taken  for  parameters  D,  C,  and  Le, 
Equations  (17)  and  (19)  reduce  to 


Tjju  +  (Pr)£  Tjj  -  -(Sc)kj,Y(AH)/C0  ,  (26) 

where  Sc  'a  the  Schmidt  number,  Sc  a^/pO,and  Pr  is  the  Prandtl  number, 
Pr  •  Sc/Ls. 

A.  Catalytic  Surface  Combustion 

In  the  absence  of  gas-phase  reaction,  k^  ■  0,  the  equations  have 
the  solutions 

Y  -  Ys  -  (Ys  -  YJ  cpm.  Sc) 

T  •  Tg  -  (Tg  -  TJ  cpOW  Pr)  (27) 


where 


.11  J\ 

cpOl,  k)  -  o(k)  Jo  exp[ -kj'fdllJdTl 


and  a(k)  is  the  constant  that  makes cp (®,  k)  »  1.  The  tabulated  values 
of  a  are  well  fitted  by 

<*(k)  -  0.570  k0’4  (28) 

The  values  Tg,  T^,  and  Y^  in  Equation  (27)  are  known;  the  value  of  Yg 
is  determined  by  the  condition  at  the  surface  T]  ■  0  as  given  in  Equation 
(21): 


from  which 


Y  «  (kg/D)  <ia/pB)1/2  Y 


Ys  -  Yw  o(Sc)/C«(Sc)  +  (kg/D)  (p/pS)1/2J 


The  heat  flux  from  the  solid  surface  by  conduction  to  the  fluid  is 


^  dy  y  m  o 


MTS  -  Tw)  a(Pr)/pB/p. 


(30) 


y.  Ln  t  a  yiizm  gyanr  tv  Tvmwqi 


Th«  .i*at  flux  to  the  catalyst  surface  by  the  exothermic,  catalyzed 
surface  reaction  is 


Qc  -  (AH)  Dd^  ■ 

S  dy  y  «  0 


(AH)  pY^/td/kg)  +  [  1/D  cr(Sc)  /£££]] 


(31) 


By  regrouping  the  transport  coefficients  we  can  rewrite  Equation  (31) 
to  read 


.1/2 


Qs  ■  o  V  /[I  +  (r/kg)  V1/2] 


(32) 


where 


r  -  D  <*(Sc)/p3o/p, 


m  -  r(AH)pY 


and  3Q  is  the  semi -empirical  constant  of  proportionality^  between  3 
and  volumetric  flow  rata  V. 

B.  Catalytic  Surface  and  Gas  Phase  Combustion 

If  reaction  in  the  gas  phase  is  included  Equations  (25)  and  (26) 
are  not  likely  to  have  an  analytical  solution.  They  are  solved  here 
numerically  for  the  temperature  and  fuel  concentration.  Functions  and 
parameter  values  used  in  the  calculations  are  given  in  Table  1.  The 
coefficients  D,  p ,  and  p,  were  taken  as  the  average  of  their  values  of 
T  -  T^  and  for  T  -  Tg. 

For  the  propane/air  system,  the  Schmidt  number  Sc  is  nearly  temper¬ 
ature  independent  with  numerical  value  unity.  The  Lewis  number  Le  is 
also  nearly  temperature  independent.  If  the  Lewis  number  is  also  taken 
to  be  unity,  Equation  (25)  can  be  replaced  by  the  following  linear  com¬ 
bination  of  Equations  (25)  and  (26)  : 


where 


•Vi +  f  ^  ‘ 0 


J  -  CT  +  (AH)Y 
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(33) 


(34) 


This  equation  has  th#  solution 

J  -  Jg  -  (Jg  -  Jm)  cpOl,  1)  (35) 

where  function  9  is  defined  for  Equation  (27). 

Values  ct  Y^  and  and  Tg  are  given  in  Equation  (21)  so  that 
may  be  computed,  while  Jg  is  known  only  when  a  value  is  assumed  for  Yg. 

In  the  overall  numerical  method  we  therefore  assume  a  value  for  the 
constant  Yg,  obtain  function  J  from  Equation  (35),  replace  Y  in  Equation 
(26)  by  its  expression  as  a  function  of  J  and  T  from  Equation  (34) ,  and 
then  solve  the  nonlinear  differential  Equation  (26)  for  T  as  a  function 
of  T).  Equation  (34)  then  gives  Y  as  a  function  of  T]  so  that  the  derivative 
Y^  can  be  computed  at  T\  ■  0  and  the  wall  boundary  condition  for  Y  in 
Equation  (21)  used  to  improve  the  assumed  value  of  After  starting 

with  the  two  extreme  values  Yg  ■  0  and  Yg  ■  Y^  we  obtained  convergence 
in  no  more  than  three  iterations  with  a  secant  method.  In  this  method, 
the  function  of  Yg  that  is  zero  according  to  the  wall  boundary  condition 
is  taken  to  vary  linearly  between  it3  two  most  recently  computed  values. 
The  resulting  linear  equation  is  solved  for  the  improved  value  of  Yg. 

The  differential  equation  to  be  solved  for  T  is  Equation  (26)  with 
Y  replaced  b>  its  expression  in  terms  of  T  and  known  function  J  as 
described  above.  The  equation  reads 

+  f  -  kG  (CT  -  J)/CP  (36) 

Multiplication  by  the  function  F(71) , 

FCT})  -  exp  [/J  fCTf  )dll/]  (37) 

gives  a  differential  equation  in  self-adjoint  form 

(F  -  G  ,  G  -  kG  F(CT  -  J)/CB  (38) 

13 

The  finite  difference  method  of  Marchuk  and  Babuska  has  fourth- 
order  accuracy,  so  that  only  a  coarse  mesh  is  needed.  It  produces  a 


l.-ll 


tri-diagonal  coefficient  matrix  from  which  a  direct  solution  can  be  ob¬ 
tained.  The  method  applies  to  linear  equations  in  self-adjoint  form 
and  not  to  Equation  (38)  because  kg  is  a  function  of  T.  We  describe  the 
method  for  the  linear  equation 

(FT^  +  ST  -  W  ,  (39) 

where  F,  S,  and  W  are  given  functions  of  T)  or  constants  independent  of  T, 
and  then  adapt  the  method  to  solve  nonlinear  Equation  (38) . 

Divide  the  interval  of  integration  Tl0  *  11  &  into  uniform  sub- 

intervals  each  of  length  h  and  let  subscript  i  denote  evaluation  at  a  point 

TL-  T1  +  ih:  F  -  F(T|,)  -  F(TL  +  ih) ,  and  similarly  for  the  other 
li  o  t  i  o 

functions.  Let  a  super  bar  denote  the  harmonic  mean: 

fl  *  6/(Fl-l/2  +  4  FIX  +  Flil/2)  ' 

Tridiagonal  matrices  B  and  M  have  for  the  three  non-zero  elements  in  their 
i-th  rows  [Boi,  Bu,  B2i]  and  [Moi,  M^]  where,  for  i  -  1,  2,..,n  on 

a  mesh  with  n  interior  points, 

>ol  ■  C2  +  0£1/2  -  0  h-l/J'** 

Bll  "  ^44  +  <Fl-l/2  '  3Fil>  Fi-l/2  +  (Fi+l/2  "  3Fi  3Fl+l/2-'/48 
»21  •  C2  +  (3F^1/2  -  F'1)  fi+1/2]/48 

Mol  "  "Fi-l/2  +  1,2  Si-l/2  “oi 
Mli  “  ^i-1/2  +  ^i+1/2  +  h  Si  Bli 

M21  ’  -f 1-1/2  +  h2  Sl+l/2  »21  <40) 

The  finite  difference  approximation  to  the  differential  equation  and 
boundary  conditions  is  then 
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where  T  ia  the  n-dimensional  vector  of  unknown  values  of  function  T  at 
the  interior  mesh  points,  W  is  the  similar  vector  of  known  values  of 
function  W,  and  £  is  the  vector  with  all  zero  entries  except  for 
(n2  Bol  W0  -  Mol  T0)  as  first  entry  and  (h2  B2n  Wfl+1  -  M2n  T^)  for 
n-th.  Since  the  matrix  M  is  tridiagonal,  symmetric  and  generally  non¬ 
singular,  Equation  (41)  uan  be  readily  solved  for  T  by  recursion.  The 

4 

truncation  erx*or  is  of  the  order  h  ,  which  allows  for  a  smaller  number 
of  mesh  points  chan  more  conventional  difference  methods.  We  tested 
this  method  over  the  interval  Os  2  with  fCH)  *  2T],  which,  for  kG  »  0 
and  appropriate  boundary  conditions,  yields  the  error  function.  The 
solution  calculated  with  n  ■  31  was  accurate  at  T)  ■  1  to  six  decimal  places. 

For  Equation  (38) ,  the  nonlinear  equation  of  physical  interest  here, 
Newton's  method  provides  a  linear  differential  equation  for  6,  the 
difference  between  the  present  guess  T°  at  the  unknown  function  T  and 
its  generally  more  correct  value  T°  +  6 . 

(F  G°  6  -  +  G°  (42) 

where  G°  is  G  evaluated  with  T  -  T°  and  G^  is  the  partial  derivative  of 
G  with  respect  to  T.  This  is  a  linear  differential  equation  for  6  in  the 
form  of  Equation  (39)  so  that  the  finite  difference  method  used  for 
that  equation  can  be  applied  here. 

A  theorem  of  Henrici^  Implies  that  the  sequence  of  improved  values 
of  T,  each  obtained  by  Newton’s  method,  will  converge  to  the  desired 
solution  if  GT  is  non-negative.  Now 

GT  "  kG,T  F(CT  "  J)/Ce  +  kGF/5 
«-kG  T  F (AH)Y/ Cf3  +  kG  F/8 

In  these  expressions  kG  is  a  positive  rate  constant  that  increases  with 
temperature,  so  that  kr  _  is  positive,  and  the  other  quantities  appearing 


ar«  all  posiclv*  so  that  the  condition  for  convergence  is  that  the  ratio 
^G,T  b®  less  than  unity.  However,  the  iterations  did  converge 

with  the  ratio  of  the  order  of  ten.  We  took  for  the  initial  guess  T°, 
the  solution  to  the  linear  equation  obtained  from  Equation  (38)  by 
setting  kQ  *  0.  The  number  of  iterations  required  was  rarely  more  than 
six  to  achieve  five  place  accuracy  in  a  solution  to  the  discrete  problem. 

To  check  the  accuracy  of:'  the  solution  we  tested  the  code  on  the 
problem 

u‘  -  2  exp(u)/(cos  l)2 
u*0atT)*»±l 

with  solution 

u  ■  2  In  [  (cos  l)/cos  T]] 

The  exact  value  of  u  at  T]»  0,  u  *-l. 231253,  was  computed  except  for  a 
6  rather  than  a  3  in  the  last  place.  The  calculation  had  31  interior 
nodes,  was  started  with  u°  ■  0,  and  converged  in  5  iterations. 


Ill  RESULTS  OF  THEORETICAL  ANALYSES 


In  che  absence  of  gas-phase  reaction  Che  results  of  the  theoretical 
analysis  show  an  interesting  relationship  between  the  exothermic  heat 
release  rate  ac  the  surface  and  the  parameters  involving  chemical  kine¬ 
tics  and  fluid  dynamics.  The  rate  of  heat  transfer  to  the  catalytic 
surface  for  a  given  fuel-air  mixture  exhibits  two  zegimea  depending  on 

the  relative  magnitude  of  the  two  terms  in  the  denominators  of  Equations 
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(31)  and  (32).  For  the  term  D  a(Sc)  (pS/p)  «  kg,  as  would  prevail 
at  low  gas  velocity  Impinging  on  the  stagnation  plate  or  with  high  rate 
constants  for  the  surface  catalyzed  reaction,  the  exothermic  reaction 
is  controlled  by  diffusive  and  convective  transport  of  reactants  to  the 
catalyst  surface.  Under  these  conditions  che  term  containing  r  is  small 
in  Equation  (32)  so  that  the  heat  release  rate  becomes  a  linear  function 
of  the  square  root  of  the  gas  velocity.  By  contrast,  at  high  convective 
flow  the  term  is  large,  and  the  heat  release  is  governed  by  surface 
reaction  kinetics  and  is  independent  of  the  gas  velocity. 

We  have  carried  out  some  calculations  for  the  C^Hg-air  system  under 
study,  using  the  model  with  surface  catalyzed  reaction  involving  Pt  as 
the  catalyst.  For  this  analysis  we  have  evaluated  the  parameters  appro¬ 
priate  to  this  reacting  system  and  to  the  physical  geometry  of  the  stag¬ 
nation  point  flow  apparatus  discussed  in  further  detail  in  the  follow¬ 
ing  section.  Using  temperature-averaged  transport  coefficients  and  the 
measured  reaction  kinetics  for  Pt-catalyzed  combustion  of  propane-air 
mixtures, ^  we  are  able  to  calculate  on  the  basis  of  Equation  (32)  the 
heat  release  rate  at  the  surface  as  a  function  of  the  propane  concen¬ 
tration  in  the  gas  stream,  the  surface  temperature  of  the  catalyst,  and 
the  reactant  flow  velocity.  The  results  shown  in  Table  2  are  limited 
to  the  surface  catalyzed  reaction  and  do  not  include  any  contribution 
from  gas-phase  combustion.  The  two  reaction  regimes  are  clearly  dis- 
tinguishabl"  by  the  use  of  parameters  in  the  table.  One  involves  the 
transport  limited  surface  process  at  low  gas  velocities,  i.e.  when  V 
«  (kg/r)  .  Under  these  conditions  the  heat  release  rate  becomes 
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proportional  to  tha  squara  root  of  tha  volumetric  flow  rata  of  tha 
reactant  gas  mixture,  i.a. ,  Qg  is  a  linear  function  of  with  slope 

a,  which  Is  nearly  independent  of  temperature.  At  higher  gas  flow  rate 
so  that  V  »  (kg/r)2,  the  heat  release  rate  is  governed  by  reaction 


When  both  gas-phase  and  surface  reaction  are  present  an  entirely 
different  pattern  emerges.  Based  on  the  results  of  the  numerical  analy 
sis  we  have  evaluated  the  relative  contribution  of  the  surface  reaction 
to  the  total  heat  release  rate  at  the  surface  (Figure  1) .  At  low  sur¬ 
face  temperature  and  high  gas  flows  the  preponderance  of  the  surface 
reaction  becomes  quite  apparent.  For  example,  with  a  gas  mixture  con¬ 
taining  1  vol  percent  CgHg  and  a  flow  rate  of  10  liters /minute,  we  find 
that  the  heat  release  traverses  from  a  predominant  surface  reaction  at 
650  K  to  a  predominant  gas-phase  reaction  at  950  K.  Reduction  of  the 
rate  constant  for  the  surface  reaction,  by  reducing  the  preexponential 
factor  by  three  orders  of  magnitude,  is  shown  to  be  accompanied  by  a 
decrease  in  the  heterogeneous  contribution  to  the  overall  heat  release 
rate  (Figure  1) . 


The  temperature  profiles  in  the  stagnation  point  boundary  layer 
offer  some  insight  into  the  effect  of  the  gas  phase  reaction  on  the 
temperature  gradient  near  the  surface.  The  temperature  profile  at  high 
gas  velocity  is  characteristic  of  the  predominantly  surface  catalyzed 
reaction  (curve  A,  Figure  2) .  At  lower  gas  velocity  (curve  B)  the 
temperature  gradient  at  the  surface  changes  from  negative  to  positive, 
as  the  gas-phase  reaction  predominates.  Similar  trends  prevail  in  the 
temperature  profiles  computed  at  the  higher  surface  temperature,  except 
that  the  gas  phase  reaction  contributes  to  the  heat  received  by  the 
surface  even  at  the  highest  reactant  flow  rate  shown.  The  temperature 
gradient  at  the  surface  is  markedly  reduced  from  that  seen  at  the  lower 
surface  temperature  for  the  same  gas  flow  rate. 

Of  further  interest  is  the  fuel  concentration  profile  in  the 
stagnation  point  boundary  layer  adjoining  the  catalytic  surface.  As 


shown  by  ehs  dsts  plotted  in  Figure  3,  eh*  prof lit  dtp ends  on  the  reletlve 
contribution  of  ges  end  eurfece  reaction.  At  the  lower  gee  velocity 
the  fuel  concentration  la  significantly  reduced  throughout  the  boundery 
layer.  This  result  is  to  be  expected,  since  the  gas  phase  contribution 
to  the  heat  release  is  much  greater  at  the  lower  velocities  (Figure  1). 
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IV  EXPERIMENTAL  DETAILS 


The  stagnation  point  flow  reactor  selected  for  our  studies  is  shown 
schematically  in  Figure  4.  A  stream  containing  the  propane-air  gas 
mixture  of  the  desired  composition  is  directed  in  an  upward  direction 
toward  the  catalytic  surface,  composed  of  a  quartz  support  plate  (10.4 
cm  x  10  cm  x  0.1  cm)  onto  which  a  series  of  12  separate  platinum  strips 
have  been  deposited  by  vacuum  evaporation  of  the  metal.  Each  individual 
strip  of  the  metal  film  (0.4  cm  x  10  cm)  extends  over  the  entire  width 
of  the  quartz  plate.  In  the  central  section  the  metal  film  over  a  dis¬ 
tance  of  8  cm  has  a  thickness  of  4  x  10  ^  cm  while  near  the  two  ends, 
over  a  distance  of  1  cm,  the  film  thickness  is  12  x  10  cm.  Each  of 
these  Pt  strips,  separated  by  a  distance  of  0.1  cm,  serves  as  a  catalytic 
surface  and  resistance  thermometer.  At  the  ends  of  each  strip  electrical 
contact  is  made  to  a  DC  power  supply  and  electrical  circuit  for  resistance 
measurement.  First  each  strip  is  heated  electrically  in  an  inert  gas 
stream  to  specified  temperatures,  as  measured  by  a  radiation  thermometer 
(with  proper  emissivity  correction,  e«  0.65),  and  its  electrical  resis¬ 
tance  recorded  at  each  temperature.  Because  of  the  difference  in  film 
thickness  the  temperature  in  the  central  portion  of  the  Pt  strip  effec¬ 
tively  governs  its  resistance.  Over  90%  of  the  length  of  each  strip  the 
temperature  was  found  to  be  uniform  within  ±  5  K  ,  with  a  sharp  temper¬ 
ature  gradient  at  the  juncture  with  the  thicker  film  ends.  The  electri¬ 
cal  power  required  to  maintain  this  temperature  is  measured  as  a  function 
of  inert  gas  velocity  impinging  on  the  plate.  Subsequently  we  determine 
the  power  required  to  keep  each  strip  at  the  same  temperate  (resistance) 
in  the  presence  of  the  fuel-air  mixture  impinging  at  the  same  total  gas 
velocity.  Because  of  exothermic  reactions  at  the  surface  and  in  the 
gas  phase,  the  difference  in  power  requirement  represents  the  energy 


£ 
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generated  by  the  chemical  reaction.  Thus,  our  meaaurementa  vere  carried 
out  at  constant  catalyst  temperatures,  but  different  gas  velocities  and 
fuel-air  mixture  ratios.  The  incoming  gas  temperature  in  this  series  of 
experiments  was  maintained  at  293  ±  5  K.  The  gas  entered  through  a  duct 
of  rectangular  cross  section  (1  cm  x  8  cm),  separated  from  the  plate  by 
a  distance  of  0.65  cm.  The  duct  was  oriented  parallel  to  the  Pt  strips. 
The  central  strip  was  located  in  the  center  of  the  nearly  two  dimen¬ 
sional  stagnation-point  flow.  During  an  initial  heating  period  in 
nitrogen  of  about  2  hours  the  film  sintered,  and  subsequently  exhibited 
little  changf  in  its  electrical  properties,  as  determined  in  a  series  of 
recalibrations  of  its  electrical  resistance  as  a  function  of  temperature. 
The  quarts/Pt  catalyst  plate  was  rigidly  suspended  inside  a  glass  cylinder 
(15  cm  in  diameter)  through  which  a  alow  stream  of  Nj  was  Paaa®d  to  sweep 
away  the  reaction  products. 


V  Eft'ERiriEuTAI.  RESULTS 


As  might  be  expected,  the  highest  rate  of  exothermic  reaction 
occurred  near  the  central  Pt  strips  exposed  to  the  stagnation  point  flow 
containing  C^Hg-air.  As  a  result  it  was  found  expedient  to  monitor 
only  the  heat  release  rates  on  the  three  central  strips,  which  were 
responsible  for  more  than  90%  of  the  heat  released  to  the  entire 
catalyst  plate  over  which  reactive  gas  was  flowing.  The  results  of  two 
experiments  carried  out  at  different  catalyst  temperatures  are  shown  in 
Figure  5  for  a  fuel-air  gas  mixture  containing  i  vol%  C^Hg.  Here  we 
correlate  the  maas.  ired  heat  release  rates  to  the  square  root  of  the 
volumetric  flow  rate  of  reactive  gas,  as  suggested  by  the  theoretical 
analysis.  Over  a  range  of  catalyst  temperatures  the  experimental  results 
are  summarized  in  Table  3  in  terms  of  the  ratio  of  the  heat  release  rate 
to  the  square  root  of  the  volumetric  flow  rate.  It  is  to  be  noted  that 
with  rising  catalyst  temperature  the  ratio  grows,  i.e.  for  a  given 
volumetric  flow  rate  the  heat  release  rate  increases.  If  we  were  dealing 
entirely  with  transport  controlled  surface  reaction  this  ratio  would 
tend  to  decline,  as  manifested  by  the  theoretical  data  listed  in  Table  3 
(column  4) .  However  the  combined  contributions  of  surface  and  gas-phase 
reactions  readily  explain  the  experimental  results  (Table  3) . 

The  theoretical  data  are  derived  from  the  model  by  setting  up  the 
heat  balance  equation  for  the  central  platinum  strip.  The  power  Qw 
required  per  unit  area  to  maintain  the  strip  at  a  fixed  temperature 
is  given  by 


Q 


w 


^A^B 


where  is  the  heat  conducted  to  the  gas  stream,  Qg  is  the  heat  lost 

to  the  surroundings  by  conduction  through  the  quartz  plate  or  in  other 
* 

ways,  and  Q_  is  the  heat  generated  by  surface  reaction.  Theoretical 

O  «  • 

expressions  are  given  for  and  Qg  ir  Equations  (30)  and  (31).  Let 
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Independent  of  flow  speed  and  temperature  and,  In  the  diffusion  controlled 
regime  prevailing,  Yg  is  always  zero.  Hence  (Jj  -  Jp^  Is  a  fixed  func¬ 
tion  of  TJ,  and  Q  is  proportional  to  dl]/dy  or  V^/2  as  shown  in  Figure  6. 

The  theoretical  calculations  of  the  separate  contributions  of  the  gas 
and  surface  reactions  to  the  total  Q  are  also  shown  in  Figure  6. 
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VI  DISCUSSION 


In  this  study  we  have  been  able  to  demonstrate  the  contribution  of 

heterogeneous  and  homogeneous  reactions  to  the  combustion  of  fuel  in  a 

stagnation  point  flow  reactor  with  a  catalytic  surface.  The  theoretical 

analysis  clearly  indicates  the  conditions  under  which  either  combustion 

mode  predominates.  At  650  K  the  surface-catalyzed  reaction  is  the  major 

component  of  the  total  heat-release  rate  (Figure  6)  over  nearly  the 

entire  range  of  fuel-air  gas  velocities.  An  entirely  different  pattern 

emerges  at  850  K.  At  this  temperature  we  note  the  gradual  transition 

from  gas-phase  predominance  to  surface-controlled  heat  release.  At  a 

3  -1 

flow  rate  of  160  cm  -s  the  heat  release  rates  by  homogeneous  and 
heterogeneous  reaction  are  equal.  At  higher  flow  rates  the  gas-phase 
contribution  declines  relative  to  the  surface  reaction.  The  general 
pattern  emerging  from  this  analysis  bears  some  similarity  to  the  nume¬ 
rical  analysis  of  CO  oxidation  catalyzed  by  Pt  in  a  monolith  reactor. 

It  is  to  be  noted  that  in  the  heterogeneous  regime  in  which  the 
heat  release  rate  is  transport  limited,  a  more  active  catalyst  will  not 
improve  the  heat  release  rate.  Only  at  impingement  gas  velocities  so 
high  that  the  catalytic  reaction  limits  the  conversion  rate  will  an 
improvement  in  catalyst  activity  be  accompanied  by  an  increase  in  the 
surface-kinetics  controlled  limit  of  the  heat-release  rate. 

The  heat  release  rate  by  catalytic  surface  reaction  attains  one 
half  of  its  maximum  value  (m  kg/r)  at  the  square  root  of  the  volumetric 
flow  rate  (kg/r).  These  values  are  given  in  Table  2.  The  flow  rates 
required  are  considerably  higher  than  employed  in  our  experimental 
studies.  Thus,  in  order  to  attain  high  combustion  efficiency  at  mode¬ 
rate  volumetric  flow  rates  of  fuel-air  it  is  necessary  to  promote  gas- 
phase  combustion. 

In  general,  the  heat  release  rates  measured  experimentally  compare 
favorably  with  the  values  derived  from  our  theoretical  model  when  the 
contribution  of  both  combustion  modes,  heterogeneous  and  homogeneous. 
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Is  taken  into  account.  The  theoretical  ratio  Q/V  1  is  independent  of 
temperature  (Table  3).  The  analytical  modeling  suggests  a  suitable 
experimental  approach  to  measure  the  contribution  of  gas  phase  reaction 
to  the  combustion  process  in  stagnation  point  flow.  The  temperature 
profiles  in  the  boundary  layer  near  the  catalytic  surface  exhibit  marked 
differences  in  the  presence  of  .homogeneous  reaction  (Figure  3) .  Thus  a 
detailed  study  of  the  temperature  distributions  at  various  volumetric 
flow  rates  of  fuel-air  and  different  catalyst  temperatures  would  allow 
elucidation  of  the  gas-phase  kinetics  and  its  contribution  in  the  overall 
combustion  process. 

In  conclusion,  the  work  described  provides  detailed  information 
on  the  contribution  of  surface-catalyzed  and  gas-phase  reactions  to  the 
combustion  of  a  fuel-air  mixture  in  stagnation  point  boundary  layer  flow. 
To  make  the  theoretical  model  more  precise  the  restriction  of  the  Lewis 
number  tc  unity  should  be  removed.  Experimentally,  temperature  probing 
within  the  boundary  layer  would  further  strengthen  our  understanding  of 
the  interplay  of  homogeneous  and  heterogeneous  reaction  kinetics  with 
the  fluid  dynamics  near  a  catalytic  surface. 
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Table  1 

FUNCTIONS  AND  PARAMETER  VALUES  FOR 
THEORETICAL  CALCULATIONS 


Quantity 

Equation 

Reference 

S0  -  3.44  x  10  3  (s  V  cm3 -rain 

(32) 

11 

AH  -  1.1  x  104  (cal-g'1) 

(8) 

17 

kg  -  1.1  x  109  exp (-17000/ RT)  (cm-s"1) 

(29) 

7 

k.  «  1.7  x  1013  exp(-21000/RT)  (cm3 -a"1- mol"1) 

(17) 

8 

p./p  *  9.36  x  10  6  T1'  (cm2*s  1) 

(15) 

* 

12 

D  -  2.70  x  10~5  T1,59  (crn^s"1) 

(16) 

16+ 

C  *  0.30  (cal-g"1^"1) 

(8) 

17 

pT  ■  0.35  (g-K*cm  3) 

(10) 

17* 

Table  3 


HEAT  RELEASE  RATES  AT 
IOUS  CATALYST  TEMPERATURES 


.  .1/2 

Q/V  cal- 

_2  _ ^ 

cm  *s  >/(cm3-s 

-V'2]  x  102 

Experimental 

Theomtlcal 

(surface) 

Theoretical 
(surface  +  gas) 

2.00 

1.72 

2.48 

2.52 

1.63 

2.47 

2.80 

1.55 

2.41 

0.98 


0.81 
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COMBUSTION  IN  A  MONOLITHIC  CATALYTIC 
COMBUSTOR  WITH  GAS  PHASE  REACTION 

by 

C.  M.  Ablow,  B.  J.  Wood,  and  K.  Wise 
SRI  International 
Menlo  Park,  California  94025 

ABSTRACT 

Experimental  observations  of  the  temperature  distributions  in  a 
multichannel  catalytic  monolith,  through  which  a  mixture  of  propane  and 
air  is  flowing,  have  been  correlated  qualitatively  with  a  theoretical 
model.  The  correlation  shows  that  gas-phase  reaction  is  initiated  by 
the  catalyzed  heterogeneous  reaction,  that  the  combustion  is  strongly 
Influenced  by  heat  conduction  in  the  web  of  the  solid,  and  that  the 
surface  reaction  rate  can  be  slow  enough  and  the  transport  fast  enough 
for  chemical  kinetic  control  of  the  process. 


I  INTRODUCTION 


The  design  of  a  multichannel  catalytic  monolith  for  use  as  a  flame 

holder^"  in  aircraft  afterburners  will  require  accurate  modeling  of  the 

initiation  of  gas-phase  combustion  by  the  heterogeneous  reaction  on  the 

catalytically  active  surface.  A  steady-state,  one-dimensional  model 

that  relates  gas  properties  average  at  each  cross  section  in  a  cylindrical 

catalyst  channel  to  the  local  wall  temperature  and  reaction  rate  by  mass 

and  heat  transfer  coefficients  has  been  formulated  with  various  degrees 

2—6 

of  approximation  by  several  authors.  Based  on  such  theoretical 
modeling,  the  predicted  behavior  of  the  system  can  be  categorized  into 
three  distinct  regimes, dependent  on  the  magnitudes  of  the  system 
variables,  such  as  gas  velocity,  temperature  and  fuel/air  ratio.  Briefly, 
these  regimes  comprise  the  following  patterns:  (1)  a  gradual  increase 
of  temperature  with  distance  into  the  monolith;  (2)  a  sharp  temperature 
increase  at  the  inlet  or  within  the  monolith;  (3)  a  major  contribution 
to  total  heat  release  by  gas-phase  reaction.  Since  the  role  played  by 
gas-phase  reactions  inside  the  monolith  is  of  special  interest,  we  have 
carried  out  a  series  of  experimental  measurements  designed  to  examine 
this  regime  in  more  detail.  Their  theoretical  interpretation  indicates 
that  axial  heat  conduction  by  the  solid  monolith  has  a  significant  effect 
on  combustor  performance,  even  in  the  central  channels  of  a  large-diameter 
monolith  where  radial  conduction  is  of  lesser  importance. 
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II  EXPERIMENTAL  STUDIES 


In  order  to  reduce  heat  losses  by  radial  conduction  and  more  closely 
approximate  the  adiabatic  conditions;  in  an  interior  channel  of  a  catalytic 
monolith,  the  single-tube  combustor  used  in  earlier  studies^  was  replaced 
by  a  multi-channel  unit.  The  need  for  an  array  of  channels  with  identical 
flows  through  them  was  clearly  established  by  our  success6  in  fitting  to 
our  model  published  data  obtained  with  a  monolith  composed  of  several 
channels^  after  taking  into  account  heat  transfer  to  the  surroundings. 

The  combustor  is  a  2.5  cm  diameter  core  cut  from  a  cordierite  mono- 

* 

lith.  It  is  14  cm  long  and  is  completely  coated  with  a  platinum- 
impregnated  alumina  washcoat.  Each  channel  has  an  open  cross  section 
0.145  cm  square  (Figure  1)  with  a  web  thickness  of  0.028  cm,  and  hydraulic 
radius  R^  -  0.0363  cm.  There  is  a  minimum  number  of  ten  channels  across 
any  diameter  of  the  core.  Three  of  the  channels  were  closed  to  gas  flow 
by  a  plug  of  Sauereisen  No.  1  cement  at  the  upstream  end  of  the  combustor, 
as  shown  in  Figure  1. 

The  catalytic  monolith  was  situated  in  a  2.5  cm  inside  diameter  fused 
quartz  tube  (Figure  2) ,  and  held  in  place  with  quartz  retaining  rings 
cemented  to  the  tube.  Propane-air  and  carbon  monoxide-air  mixtures  were 
prepared  by  metering  the  gases  through  flow  meters.  The  selected  gas 
mixture  entered  the  fused  quartz  tube  and  flowed  through  a  30  cm  long 
section  heated  by  an  electri^  furnace.  At  the  end  of  this  preheat  section 
the  combustible  gas  mixture  entered  the  combustor,  which  was  also  heated 
by  an  electric  furnace.  At  the  outlet  of  the  catalytic  monolith,  the 
gases  passed  through  a  15  cm  long  section  of  the  fused  quartz  tube  before 
being  vented  into  a  fume  hood. 

The  axial  and  radial  steady-state  wall  temperature  profiles  in  the 
closed  channels  of  the  combuster  were  monitored  by  a  set  of  three  chromel- 
alumel  thermocouples  in  0.1  cm  sheaths  with  exposed  junctions.  A  fourth 
thermocouple  was  located  in  a  central  open  channel  to  monitor  the  gas 
temperature.  These  thermocouples  were  attached  to  a  single  movable 


*The  material  was  kindly  provided  by  Dr.  L.  L.  Hegedus,  General  Motors 
Research  Laboratories,  Warren,  Michigan. 


fixture  so  that  they  could  be  moved  axially  as  a  unit  to  any  desired 
cross  section  within  the  combustor.  An  additional  thermocouple  was  fixed 
in  the  gas  stream  1.5  cm  ahead  of  the  catalytic  monolith. 

Samples  of  gas  from  the  feed  stream  and  the  product  stream  of  the 
reactor  were  taken  with  gas-tight  syringes.  The  samples  were  analyzed 
for  propane,  carbon  monoxide,  and  carbon  dioxide  by  injecting  them  into 
a  gas  chromatograph  with  Porapak  T  and  5A  Molecular  Sieve  columns.  Com¬ 
plete  conversion  of  the  fuel  to  CO,  occurred  during  passage  of  the  gas 
through  the  catalyst  in  al  runs  with  V  <  8  i-min  .  With  V  ■  8  i-min  , 

96  vol%  of  the  propane  in  the  feed  stream  was  found  to  be  consumed. 

Results  of  the  experiments  are  presented  in  Figures  3  through  5. 

The  radial  profiles  of  the  temperature  at  various  axial  distances  from 
the  inlet  to  the  monolith  (Figure  4)  exhibit  only  small  gradients.  Hence, 
neglect  of  radial  heat  transfer  in  our  analytical  model  represents  a 
reasonable  approximation  of  the  exper 'mental  situation. 


Ill  THEORETICAL  MODEL 


* 

A.  Duct  Flow  Model 

The  steady-state  temperature  and  concentration  distributions  in  each 
duct  cross  section  are  modeled  by  their  average  plug  flow  values  in  the 
gas  phase,  denoted  by  subscript  G,  and  their  values  at  the  catalyst  sur¬ 
face,  subscript  S.  The  heat  balance  in  a  cross  section  of  the  fluid 
reads 

d2T 

^  dx  (PGV  CG  V  ^G  ~2 

dx 

hT^TG  ”  V  +  RH^k  GPG  YG  ^ 

where,  for  the  fuel-lean  case,  the  nearly  constant  concentration  of 
oxidizer  is  included  in  the  reaction  rate  constant  k^.  The  heat  transfer 

yj 

coefficient  may  be  expressed  in  terms  of  either  the  Stanton  or  Nusselt 
numbers : 


hT  *  PGV  CGSt  "  (V4RH)NU 

The  fuel  mass  balance  in  the  fluid  gives 

j  d2(prY  ) 

i  '  *h D  — t— 

dx 

-  Vc  *0  -  psV  -  Wo  yg  •  <2) 

The  mass  transfer  coefficient  h^  may  be  written  in  terms  of  the  diffusion 
coefficient  D  of  the  fuel  through  the  gas  mixture  and  the  Sherwood  number 
Sh: 


h^  -  (D/4RH)Sh 


The  notation  used  in  this  paper  is  identical  to  that 
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found  in  reference  6. 


The  fuel  Chat  diffuses  Co  the  surface  of  Che  caCalysC  is  consumed 
by  Che  reacClon: 


hM(pGYG  ~  PsV  "  kSPSYS 


(3) 


Finally,  heac  balance  In  Che  ducc  wall  may  be  wriccen  as: 


d3Te 

^w^w  ,  2 
dx 


+  hi<T0  -  V 


+  QVsYs  ■  »„  VTs  -  V 


w 


where  N  Is  Che  racio  of  Che  heac  cransfer  coefflcienCs  for  exterior  and 
w 

interior  temperature  differences  (subscripts  E  and  S,  respectively). 


B.  Governing  Parameters 

The  governing  parameters  are  found  by  a  nondimens ionalizat ion  of 
Che  equations.  The  dimensionless  distance  and  temperature,  §  and  T,  are 
taken  to  be 

§  -  x  St/Ry  ,  T  -  T/ATa  (5) 

where  the  adiabatic  temperature  rise  AT^,  is  defined  by 

4TA  •  0  WCG 

*v 

and  a  constant,  uniform,  average  value  is  taken  for  the  Stanton  number, 
St,  and  for  the  other  transport  parameters,  Hu,  Sh,  D,  \  q,  etc. 

Since  m  *  p^v  is  a  constant,  Equation  (1)  may  be  divided  by  mCg(St)ATA 
to  obtain 


T'o  ■  rV*%  -  (tg  -  V 
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where  Pe^  Che  Peclet  number  for  heat  transfer  in  Che  scream.  Is  defined 
by 


Peh  ‘  *B  W  WC 

and  Che  primes  denote  differentiation  with  respect  to  §. 

Introduction  of  the  unbumed  fuel  fraction  y  as  a  dependent  variable, 
and  use  of  the  constant  pressure  relation,  p„T,,  -  constant,  reduce  the 

u  U 

equation  to 

Tg'  "  TG^Peh  ”  <TG  "  TS^  +  *G  yG^TG  * 

^G  *  kG  PGOTGO//°^St^TA  ^ 

Similar  manipulation  of  Equation  (2)  results  in 

*0  "  yG  /Pem  “  (Sh/Nu  U)CyG  *  (pS/pG)yS]  ’  Vg/tG 

where  Pe  ,  the  Peclet  number  for  mass  transfer,  is  defined  by 
m 

Pem  "  *H  m*GD  SE  *  Lepeh 

and  Le,  the  Lewis  number,  is  defined  by  Le  -  X/p^C^D.  The  analogy  bet¬ 
ween  heat  and  mass  transfer,  Sh  ■  Nu,  gives 

y*G"  ‘yG/Pem  '  &G  ’  (p3/pG)yS]/(Le)  "  Vg/tG  *  (7) 

Equation  (3)  reduces  by  the  same  method  to 

yG  -  (pg/pG)ys  ■  <L«)^sys/,Ts  ’ 

"  kSPG0TG0^St^TA  ^ 
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Finally  Equation  (4)  becomes 


Ts/P'»  *  tg  -  Ts 


+  VS/TS 


VTS  -  te) 


where  Pe^,  the  Peclet  number  giving  the  ratio  of  heat  transfer  between 
gas  and  solid  to  conduction  in  the  solid  wall,  is  defined  by 

■  rh2  yv»<st)2  • 

The  fuel  fraction  yg  at  the  catalytic  surface  can  be  eliminated  from 

Equations  (6),  (7)  and  (9)  by  use  of  Equation  (8).  There  remain  the  three 

following  equations  for  the  three  variables  y_,  t„,  and  T„: 

b  b  S 


T'o/Peh  +  T'  G  ■  -(TG  -  TS>  +  kG  yG/TG 


-  +  K  '  -»e  -  Vg/tg 


TS/Pew  +  TG  '  TS  +  %  ‘  VTS  '  V  ' 


F  -  1/ (Le  +  TG/ks) 


fraction  F  may  be  recognized  as  the  factor  that  indicates  the  control  of 
the  heterogeneous,  catalyzed  reaction  either  by  chemical  kinetics,  if  kg 
is  small  so  that  the  second  term  dominates,  or  by  diffusion,  if  kc  is 

i  8 

large. 

This  system  of  differential  equations  has  a  definite  solution  only 
if  six  I. -•••it.dary  cur.  ’  .ions  are  specified,  such  as  the  following.  Entry 
values  of  gas  composition  and  temperature  determine  y^  and  at  j  ■  0. 
Little  heat  is  conducted  through  the  ends  of  monolith  duct  walls  so  that 
one  may  take 

T1,,  ■  0  at  J  »  0  and  §  ■  L 
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The  fuel  concentration  and  temperature  in  the  gaa  stream  at  §  ■  L,  the 
exit  from  the  monolith,  will  be  continuous  with  conditions  downstream. 
Precise  conformance  to  these  boundary  conditions  is  not  needed  to  apply 
the  model  qualitatively. 

C.  Application  of  Theoretical  Model 

Numerical  values  of  the  parameters  were  used  to  further  simplify 

9 

the  model  in  the  present  application.  Theoretical  calculations  for  the 
viscosity,  conductivity,  and  specific  heat  of  the  pure  gases  ^  and  Oj 
have  been  averaged  to  give  the  values  presented  in  Table  1.  The  volume¬ 
tric  flow  rate  V,  an  experimental  variable,  together  with  the  density 
from  the  perfect  gas  law,  determine  the  Reynolds  number.  The  diffusivity 
of  propane  has  been  taken  to  be  the  same  as  that  for  carbon  dioxide, 
a  gas  of  the  same  molecular  weight .  The  calculated  value  of  the  Lewis 
number  used  in  the  model  is  near  unity  at  all  temperatures. 

*  _1 

Experimental  values  of  the  flow  rate  V  were  no  greater  than  10  f-min  . 

Hence,  the  Reynolds  number  is  small  and  a  laminar  flow  value  for  the 
Nusselt  number^  can  be  used  to  compute  the  Stanton  numbers  tabulated 
in  Table  1. 

The  stream  Peclet  numbers  for  mass  and  heat  are  equal  in  a  gas  with 
Le  »  1  and  may  be  readily  computed, as  can  the  wall  Peclet  number  when 
the  wall  conductivity  is  known.  Our  experimental  determination  of  the  wall 
conductivity  is  described  in  the  Appendix.  Since  the  wall  Peclet  numbers 
are  two  orders  of  magnitude  smaller  than  the  stream  numbers  and  since 
these  numbers  appear  as  reciprocals,  there  is  some  justification  for 
simplifying  the  equations  of  the  model  terms  containing  Pem  and  Pe^. 

Finally  for  a  duct  in  the  center  of  the  monolith,  surrounded  by 
ducts  containing  similar  reactive  flows,  radial  heat  transfer  can  be 
neglected.  Then  Nw  ■  0  and  the  equations  of  the  model  reduce  to: 

TG  "  “  (TG  -  TS)  +  \  yG/TG  (13) 

y*G  "  ‘  FyG  *  \  yG/TG  (14) 
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TS/Pew  +  TG  "  TS  +  "  ° 


(15) 


A  linear  combination  of  the  three  equations, 

TS/Pew  -  tG  '  y'o  *  0 
can  be  integrated  to  give 

TS,Pev  -  TG  +  yG  '  TG0  -  1  (16) 

Where  the  constant  of  Integration  has  been  evaluated  by  use  of  the 
boundary  conditions  at  §■  0.  Differential  equations  (13),  (14),  and  (16) 
describe  the  simplified  model.  Boundary  conditions  for  this  third-order 
system  may  be  taken  to  be 

yG  '  1  >  TG  '  TGO  aC  5  *  0 
r's  *  0  at  5  *  L 

Since  second  derivative  terms  in  yG  and  have  been  dropped,  no  down¬ 
stream  boundary  conditions  are  needed  for  these  variables. 
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IV  DISCUSSION 


A.  Effects  of  Flow  Velocity  on  Temperature  Distribution 

The  steady  state  wall-temperature  profiles  observed  in  the  monolith 
at  various  flc  ra;es  (Figure  3)  demonstrate  a  pronounced  variation  in 
shape  as  the  inlet  gas  velocity  is  changed.  Since  the  dimensionless 
distance  variable  §  in  the  model  (Equation  5)  is  proportional  to  V,  one 
expects  temperature  distributions  at  different  speeds  that  are  similar 
to  one  another,  differing  only  in  the  distance  scale.  Some  evidence  of 
this  effect  is  the  experimental  observation  that  the  temperature  maxima 
on  the  curves  move  away  from  the  entry  as  the  flow  velocity  is  increased. 
However,  the  distance  scale  change  does  not  explain  the  difference  in 
the  values  of  the  temperature  maxima  or  other  differeices  in  the  shapes 
of  the  curves.  These  effects  can  only  appear  in  the  model  through  the 
Peclet  parameters  which  are  strongly  velocity  dependent.  Thus  heat  con¬ 
duction  has  an  important  effect  on  combustion.  As  the  flow  velocity 
tends  to  2ero,  so  do  Pew  and  the  second  derivative  of  Tg  [Equation  (5)]. 
The  model  therefore  predicts  that  the  variation  in  the  surface  temperature 
with  distance  will  become  more  nearly  linear  as  the  velocity  decreases. 

The  conclusion  is  supported  by  the  data  in  Figure  3. 

B.  Contribution  of  Surface  Reaction 

Let  J  be  an  average  value  for  (F  +  k^/T^)  a  short  distance 

from  the  inlet  of  the  monolith.  Then  Equation  (14)  can  be  approximated 
by 


y 


G 


The  solution  of  this  equation  that  satisfies  the  initial  condition  yQ  -  1 
at  §  ■  0  is 


y 


G 


% 


2-11 


Data  in  Figure  3  show  that  over  the  range  of  V  employed, the  surface 
temperature  reaches  a  maximum  some  distance  from  the  inlet  to  the  mono¬ 
lith,  indicative  of  a  significant  unbumt  fuel  fraction  y^  in  that  region. 
Baaed  on  the  model,  this  rise  in  temperature  within  the  monolith  can 
occur  only  if  ?  <  1,  where  is  the  value  of  §  for  x  »  1.  Now 

h  '  St/RH 

*  2/(0.0363  V) 

based  on  the  parameter  values  given  in  Table  1.  Since  V  £  8  A /min ,  this 

gives  s  6  and  F  <  1/6.  Such  small  values  of  *F  imply  kinetically 

controlled  conditions  at  the  surface  with  kg  of  the  order  of  Tg/6. 

Apparently  the  catalytic  surface  activity  is  about  three  to  four  orders 

13 

of  magnitude  less  than  that  measured  for  Pt  wire. 

C.  Contribution  of  Gas  Phase  Reaction 

The  temperature  difference  (Tg  -  TQ>  plotted  in  Figure  5  is  propor¬ 
tional  to  the  dimensionless  variable  A  «  Te  -  T„.  Elimination  of  T_  from 

b  U  (j 

the  transport  terms  of  Equation  (10)  gives 

(A"  -  Tg)/Pe^  +  Tg  -  A'  -  A  »  kG  Yq/t g 

At  the  maximum  value  for  the  surface  temperature,  Tg  »  0  and  Tg  is 
negative.  The  experimental  results  (cf.  Figures  3  and  5)  demonstrate 
also  that  A  is  near  zero  or  negative,  A'  is  negative  and  A"  positive,  so 
that  each  term  on  the  left  of  the  equation  is  positive  or  zero.  A 
positive  gas  reaction  term  is  needed  to  balance  the  equation.  Thus  the 
model  and  the  data  correlate  only  if  gas  phase  reaction  is  present  at 
least  in  the  vicinity  of  the  temperature  maxima  in  the  catalyst. 

If  thermal  conduction  is  neglected  in  both  the  gas  and  the  wall, 
the  model  has  an  analytic  solution.  The  solution  exhibits  for  Le  a  1  a 
surface  temperature  that  increases  monotonically  with  distance  into  the 
monolith,  a  condition  irreconcilable  with  our  experimental  data. 
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In  the  following  section  (Technical  Report  No.  3)  we  apply  numerical 
analysis  to  the  model  and  compare  the  theoretical  data  with  the  experi¬ 
mental  results. 


! 

i 
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V  CONCLUSIONS 


Our  experimental  multichannel  catalytic  combustor  presents 
sufficiently  well-defined  conditions  for  theoretical  modeling.  Temper¬ 
ature  profiles  measured  In  the  channels  of  the  monolith  correlate 
qualitatively  with  a  model  that  distinguishes  between  temperatures  or 
fuel  fractions  in  the  gas  and  at  the  catalytic  surface. 

The  correlation  indicates  that  axial  heat  conduction  in  the  monolith 
walls  has  an  important  Influence  on  the  pattern  of  reaction  in  the  ducts 
of  the  monolith,  that  the  heterogeneous  reaction  on  the  catalytic  sur¬ 
faces  of  the  monolith  may  be  kinetlcally  controlled,  and  that  the 
catalyzed  reaction  of  the  duct  walls  can  act  to  initiate  gas-phase 
reaction. 
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FIGURE  1  END  VIEW  OF  MONOLITHIC  CATALYTIC  COMBUSTOR 
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FIGURE  2  SCHEMATIC  DIAGRAM  OF  EXPERIMENTAL  APPARATUS 
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DEPENDENCE  OF  COMBUSTOR  WALL  TEMPERATURE 
ON  AXIAL  DISTANCE  AND  GAS  VELOCITY 

(Inlet  Conditioni:  Propane/Air  Ratio:  0.006;  Gal  Temp: 
220  ±  10°C) 
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Appendix 

THERMAL  CONDUCTIVITY  OF  THE  MONOLITH 


If  the  temperature  of  the  surface  of  a  semi-infinite  slab  is  raised 
at  a  rate  proportional  to  tn/^,  where  t  is  the  time  and  n  is  an  integer, 
then  the  temperature  rise  within  the  slab  of  depth  x  is  proportional 

to  in  erfc  (x/Vxt)  where  x  is  the  thermal  diffusivity  in  the  slab  and 

n  lc 

i  erfc  is  the  n-th  integral  of  the  complement  of  the  error  function.  J 

One  may  show  that  the  dominant  term  of  an  expansion  of  In  [ in  erfc  (u)] , 

2 

is  (-u  )  for  any  n  as  u  becomes  large. 

To  determine  x,  the  thermal  diffusivity,  for  the  monolith  core  used 
in  our  experiments,  the  furnaces  (Figure  2)  were  repositioned  so  that 
the  juncture  between  them  coincided  with  the  location  of  the  inlet  end 
of  the  monolith.  With  the  system  at  room  temperature  (298  K)  *and  without 
any  gas  flow,  the  rreheat  section  furnace  (Figure  2)  was  energized, 
causing  the  temperature  of  the  inlet  face  of  the  monolith  to  rise  mono- 
tonocially  with  time.  The  face  temperature,  T^,  was  monitored  with  a 
thermocouple  situated  at  the  mouth  of  the  central  blocked  channel  of  the 
monolith.  A  second  thermocouple,  located  in  the  adjacent  channel  exactly 
3.8  cm  from  the  inlet  face,  measured  T2  simultaneously.  The  temperatures 
recorded  during  the  initial  eight-minute  interval  are  tabulated  in 
Table  Al. 

A  plot  of  ln(T^  -  298)  versus  In  T  shows  that  the  front  face  temper¬ 
ature  rose  at  a  rate  proportional  to  tn  with  n  *  4.5.  Although  n  is  not 
an  integer,  it  may  be  expected  that  In  (T2  -  298)  will  be  a  linear 
function  of  t  ^  for  small  t.  The  data  from  t  •  5  to  8  minutes  does 
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follow  such  a  linear  relation.  From  the  slope  one  finds  x  *  0.11  cm  *8 

The  average  bulk  density  of  the  monolith,  evaluated  from  the  mea- 

-3 

sured  mass  and  dimensions  of  a  cylindrical  section  is  p  *  0.44  (g-cm  ). 
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This  value  is  in  reasonable  agreement  with  published  data  for  the 
solid  density  of  refractory  glass  ceramic  material  based  on  cordierite, 


after  correcting  for  a  volumetric  void  ratio  of  0.70.  The  specific  heat 
for  this  material  is  Cw  ■  0.21  (cal*g  Based  on  these  parameters 

we  find 


Technical  Report  Ho.  3 
NUMERICAL  ANALYSIS  OF  CATALYTIC  COMBUSTION 

by 

C.  M.  Ablow  and  H.  Vise 

ABSTRACT 

A  numerical  analysis  has  been  carried  out  of  the  contribution  of 
surface-catalyzed  and  gas-phase  reactions  to  the  combustion  of  a  dilute 
propane-air  mixture  in  a  catalytic  monolith  reactor  containing  highly 
dispersed  platinum  on  its  walls.  The  results  obtained  are  in  satisfac¬ 
tory  agreement  with  the  experimental  measurements  of  the  temperature 
profiles  in  the  central  duct  of  a  multichannel  monolith  reactor. 
Specifically  the  numerical  investigation  demonstrates  the  effects  of 
inlet  gas  velocity  and  inlet  surface  temperature  on  the  location  of  the 
temperature  maximum  in  the  reactor  and  the  monolith  length  required  to 
attain  a  specified  degree  of  conversion  of  fuel  in  the  gas  stream. 


I. 


Introduction 


Catalytic  combustion  in  a  tubular  reactor  has  been  modeled  by  several 
authors.1'^  The  one-dimensional  model  with  slug  flow  appears  to  provide 
a  suitable  approximation  when  account  is  taken  of : 

a  Reaction  kinetics  on  the  catalytic  wall  and  in  the  gas 
stream. 

e  Convection  of  heat  and  mass  by  the  stream 

a  Axial  conduction  of  heat  in  the  stream  and  in  the  tube 
wall 

e  Axial  diffusion  of  reactants  in  the  stream 

a  Radial  diffusion  of  reactants  and  conduction  of  heat  in 
the  stream 

a  Radial  conduction  of  heat  through  the  tube  wall  to  the 
exterior. 

In  order  to  obtain  a  better  understanding  of  the  interplay  of  the 
fluid  dynamic  and  kinetic  parameters  and  their  relative  importance  in 
monolith  combustion,  we  have  carried  out  an  analytic  and  numerical  in¬ 
vestigation  of  the  solutions  to  a  model  applicable  to  the  experimental 
studies  presented  in  the  preceding  Technical  Report  No.  2.  Again  the 
contributions  of  surface-catalyzed  and  gas-phase  reactions  are  of  special 
interest  in  our  analysis  of  monolith  combustion. 


II.  The  Model 

For  a  specified  gas  mixture  of  fuel  and  oxidizer  the  adiabatic 
temperature  rise  AT  is  a  suitable  unit  for  normallizing  Che  temperature. 

A 

Nondimens ional  temperatures  of  the  gas  and  the  catalytic  surface,  Tg  and 
Tg ,  are  therefore  defined  by: 


tg  ■  V6IA  •  Ts  ■  V4IA 


where  TQ  and  Tg  are  the  actual  temperatures.  The  ratio  yQ  or  yg  of  the 
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fuel  concentration  Y^  or  Yg  to  the  initial  concentration  in  the  gee 
stream  Y^q  is  a  nondlmenslonal  fuel  concentration  variable,  so  that 

yG  "  YG/YGO  *  y3  "  YS/YGO 
The  nondlmenslonal  axial  distance  £  is  taken  to  be 

l  -  X  St/Ry 

where  x  Is  the  actual  distance  from  the  monolith  entry,  is  the  hydrau¬ 
lic  radius  of  the  duct,  and  St  is  the  Stanton  number,  the  ratio  of  heat 
convection  to  heat  transfer  between  wall  and  stream.  If  follows  that 
a  temperature  difference  between  two  points,  separated  by  one  nondimen- 
sional  unit  on  the  stream  axis,  will  produce  the  same  heat  flux  as  that 
difference  would  generate  between  two  points  in  the  same  duct  cross- 
section,  one  on  the  axis  and  one  at  the  wall. 

The  steady  state  equations  for  heat  and  mass  balance  in  the  gas 
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stream  and  at  the  duct  wall  can  be  combined  into  the  three  following 
nondlmenslonal  differential  equations: 

TG/Peb  ‘  +'c  "  (TG  -  TS>  +  kG  yG/TG  "  °  *  (1) 

yG/LePeh  -  yG  “  FyG  -  kG  yG/TG  "  0  ’  (2) 

TS/PeW  +  TG  ™  TS  +  ^G  *  \  <TS  ’  TE)  ’  (3) 

F  -  l/(Le  +  TG/ks) 

where  the  primes  denote  differentiation  with  respect  to  5 ,  Le  is  the 
Lewis  number  (the  ratio  of  heat  conduction  to  fuel  diffusion) ,  Pe^  the 
Pec.let  number  (the  ratio  of  thermal  convection  to  conduction  in  the  gas 
stream-! ,  and  Pew  the  mixed  Peclet  number  (the  ratio  of  thermal  convection 
in  the  stream  to  conduction  in  the  duct  wall)  ,  constant  is  the  ratio 
of  the  heat  transfer  coefficients  from  the  wall  to  the  exterior. 
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The  nondimens ional  kinetic  rate  constants  for  gas  phase  reaction  and 
surface  reaction  kg  are  taken  to  be  constant  multiples  of  Arrhenius 
exponential  functions  of  the  temperature,  an  adequate  approximation  for 
fuel  lean  conditions  (cf.  Technical  Report  No.  2). 

Because  the  monolith  is  of  finite  size  the  duct  wall  conducts  no 
heat  through  X  *  0  and  X  *  L  so  that 

Tg  *  0  at  5  ■  0  and  5  *  '  W 

For  a  duct  in  the  center  of  the  multichannel  monolith,  radial  heat 
transfer  to  the  exterior  is  negligible,  *  0.  Continuity  with  in¬ 
coming  and  outgoing  unconfined  streams,  in  which  radial  heat  transfer  is 
also  negligible,  makes  a  reasonably  accurate  model  of  the  monolith  envi¬ 
ronment.  As  a  result  one  obtains 


TVPS  ■  t'g  +  h  Vtg  *  0 

(5) 

!'g/u  Peh  -  *e  -  ‘‘o  Vtg  -  0 

(6) 

valid  for  %  <  0  or  §  >  together  with  the  boundary  conditions 

TG  ■  TGO  •  PG  •  1  at  5  "  ’ 
t' '  ■  y _  -  0  at  §  »  ® 

0  Cj 

tG»  t'g,  yG,  and  yG  are  continuous  at 

5*0  and  §  *  (7) 

In  the  next  section  we  make  further  approximations  and  simplifications 
of  the  model  by  consideration  of  the  relative  magnitudes  of  the  transport 
parameters. 
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Ill  Application 


In  the  experimental  set-up  of  monolith  combustion  described  in 
Technical  Report  No.  2,  the  configuration  of  the  monolith  was  so  designed 


as  to  reduce  radial  heat  losses  to  a  minimum  for  the 

ducts,  as  confirmed  by  the  temperature  measurements 
the  monolith.  With  *  0,  the  sum  of  Equations  (1) 

centrally  located 

in  a  cross  section  of 

,  (2) ,  and  (3)  reads 

(TG  +  yG/Le)/Peh  +  TS/PeW  "  TG  +  * 

This  may  be  integrated  to  give 

<tG  +  VLe)/Peh  +  TS/PeW  "  TG  +  yG  "  C 

(4) 

where  c  is  a  constant  of  integration.  Based  on  the  transport  parameters 
listed  in  Table  1  one  finds  that  Pe^  differs  by  two  orders  of  magnitude 
from  Pe^.  It  is  therefore  reasonable  to  neglect  the  terms  containing 

Pe^  in  Equation  (4) .  Also  wiuh  further  approximation  we  neglect  these 
terms  in  Equations  (1)  and  (2)  to  obtain  the  first  order  system 

T'g  ‘  TS  ■  TG  +  *0  VT0 

(8) 

y’a  ■  -<F  +  VTG,yG 

(9) 

TS  '  PaW  (tG  +yG  ‘  =  > 

(10) 

The  neglect  of  thermal  and  mass  diffusion  in  the  stream  allows  dis¬ 
continuities  of  the  diffusive  flux  at  §  *  0  and  5  ■  5^  so  that  the 
boundary  conditions  can  be  simplified  to 

TG  ■  T00  •  *8  ■  1  •  T's  -  0  at  5  -  0 

TS  '  0  «  5  - 

The  conditions  at  §  ■  0  show  that 
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r-3SraSS3tBesatsra 


c  ■  1  +  T 


GO 

Equations  (8),  (9),  and  (10)  have  been  solved  numerically  by  the  use 
of  a  standard  computer  routine  for  solving  initial -value  problems.  To 
use  the  routine  one  specifies  a  value  for  T,  at  §  *  0.  The  correct  rc 
has  been  chosen  if  the  unused  boundary  conditions,  t's  *  0  at  ?  «  is 
found  to  be  satisfied.  However,  the  numerical  solutions  show  the  surface 
temperature  increasing  to  a  maximum  and  then  falling  with  no  leveling 
off  to  meet  the  downstream  conditions.  This  behavior  is  probably  true 
for  linearized  versions  of  the  equations  which  have  solutions  for  rc  of 
the  form 


>:• 


l 

\ 

i 

I' 


f 


A  +  Be"b5  +  Ce"c5  + 


De 


-d£ 


where  A,  B,  C,  D,  b,  c,  and  d  are  real  constants.  The  equation  t"  ■  0 

9 

may  then  be  written 


bB  +  cC  e 


(b-c)| 


+  dDe 


(b-d)5 


This  equation,  like  a  three-term  polynomial,  can  have  at  most  two  roots. 

Thus  T*  can  be  zero  in  at  most  two  points.  One  of  these  has  been  fixed 
o 

at  §  ■  0.  The  other  is  at  the  maximum  for  Tg.  The  downstream  boundary 
condition  cannot  be  satisfied  in  the  linearized  cases  and  probably  not 
for  che  nonlinear  equations. 

Using  the  coefficients  listed  in  Table  1,  we  have  obtained  numerical 
solutions  for  the  axial  temperature  distribution  in  the  catalytic  duct  wall 
as  a  function  of  the  inlet  surface  temperature,  TgQ  (Figure  1).  It  will 
be  noted  that  at  the  fixed  inlet  gas  temperature, the  value  of  TgQ  has  a 
pronounced  effect  on  the  temperature  profile.  Also  a  change  in  the 
surface  reaction  rate,  brought  about  mathematically  by  reducing  the 
preexponential  coefficient  by  several  orders  of  magnitude  brings  about 
a  marked  change  in  the  location  of  the  temperature  maximum  and  its  ab¬ 
solute  value  (Figure  2) .  At  higher  inlet  gas  velocities  (V  ■  6  1/min) 
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Che  value  of  Che  inlec  surface  CemperaCure  has  even  a  more  pronounced 

effect  on  the  temperature  profile  (Figure  3).  At  an  inlet  temperature 

Tso  slightly  above  the  inlet  gas  temperature  (TQ0  -  TSQ  *•  1.5  K)  the 

surface  reaction  rate  is  Coo  slow  to  cause  significant  heat  release.  At 

somewhat  higher  inlet  surface  temperatures  a  marked  exothermic  reaction 

is  initiated.  Under  these  conditions  we  observe  that  the  maximum  in  the 

temperature  profile  moves  toward  the  duct  inlet  as  is  increased. 

This  effect  of  the  inlet  surface  temperature  on  the  location  of  the 

temperature  maximum  X  i3  further  demonstrated  in  Figure  4. 

max 

Of  special  interest  is  che  double-value  property  of  X  ,  solutions 
in  which  either  the  wall  -  temperature  is  very  near  the  gas  temperature  at 
entry  and  little  heat  is  released,  on  one  in  which  the  two  temperatures 
are  quite  different  and  a  great  deal  of  reaction  has  occurred.  For 
example  at  V  ■  6  1/min  a  maximum  in  the  temperature  profile  within  the 
catalyst  duct  can  be  established  at  a  distance  of  3  cm  from  the  inlet 
to  the  duct  when  T  _  ■  504  K  and  T0rt  *  539  K.  However  the  maximum 
temperatures  attained  at  the  duct  wall  correspond  to  505  K  and  753  K, 
respectively,  as  can  be  seen  from  the  curves  in  Figure  5,  in  which  we 
have  correlated  the  maximum  surface  temperatures  in  the  duct  with  the 
inlet  surface  temperatures.  These  two  steady-states  may  bs  the  result 
of  two  reaction  regimes,  one  a  predominantly  surface-catalyzed  reaction 
with  limited  heat  release  rate,  and  the  other  a  surface  initiated  gas- 
phase  reaction  with  high  heat  release  rate. 

As  expected  this  increase  in  heat  release  rate  is  reflected  in  the 
fuel  composition  profiles  obtained  at  several  inlet  surface  temperatures 
(Figure  6).  Again  we  note  the  significant  reduction  in  the  fraction  of 
residual  fuel  as  we  increase  the  value  of  Tg^. 

We  have  attempted  to  correlate  the  numerical  results  with  the 
experimental  data  of  the  preceding  Technical  Report  No.  2.  In  Figures  2 
and  3  we  show  the  experimental  temperatures  measured  in  the  catalytic 
monolith  at  two  inlet  gas  velocities.  In  fitting  the  experimental  data 
to  the  theoretical  curves  we  obtain  a  "best  fit"  on  reducing  the  pre¬ 
exponential  factor  in  the  Arrhenius  equation  for  the  surface  catalyzed 
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reaction  by  a  factor  of  10  4  from  that  given  in  Table  1.  The  value  of 
k  listed  was  obtained  in  separate  experiments  employing  a  Pt-ribbon 
l-ather  than  the  type  of  wall  coating  with  dispersed  Pt  on  an  insulator 
support  employed  in  the  current  study.  Most  likely  the  catalyst  in  the 
monolith  had  lost  some  of  the  activity  because  of  surface  contaminants 
introduced  during  its  preparation  or  during  the  course  of  experimentation 
on  exposure  to  feed  gas  (propane  and  air)  of  modes  or  unknown  purity. 

No  attempt  was  made  to  "clean"  the  catalytic  surface  preceding  an  ex- 
periment.  It  is  of  interest  that  little  change  in  activity  was  detect¬ 
able  over  the  course  of  this  experimental  study. 

Using  the  reduced  value  for  catalyst  activity,  we  observe  a  satis¬ 
factory  correlation  between  experimental  results  and  numerical  data. 

We  conclude  that  the  gas-phase  reaction  initiated  by  the  surface- 
catalyzed  reaction  makes  a  significant  contribution  to  the  overall  com¬ 
bustion  process.  By  this  mode  of  operation  in  a  catalytic  duct  the  l.eat 
release  rate  can  be  greatly  augmented  and  the  reactor  volume  considerably 
reduced.  As  a  result  of  the  homogeneous  and  heterogeneous  reactions, 
the  high  heat  release  rates  can  provide  a  suitable  source  for  flame 
holding  in  an  afterburner. 
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CONSTANTS  AND  FUNCTION  COEFFICIENTS 
AVERAGE  ABOUT  T  -  600  K 

(for  the  experimental  conditions  of 
Figure  1) 


AT.  «  358  K 
A 

St  »  2.18/V 
R^j  *  0.0363  cm 
Le  -  1.0 

-  0.225  V2 
Pew  -  0.0031  V2 

£s  -  8.67  x  107  exp(-8500/T<,) 

-  2.15  x  105  exp (-10500 /Tg) 
tgq  -  500  K 

Nw-° 

L  ■  15  cm 

'  _i 

In  the  above  V  has  the  units  (l*rain  )  and  T(K) . 
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FIGURE  1 


TEMPERATURE  PROFILE  IN  CATALYTIC  MONOLITH  DUC” 

(Inlet:  gat  temperature  »  500  K,  0.5  vol.  %  C3Hg  in  air;  surface  rate  constant  k" 
8.67  x  104  exp  (-8500/T.);  k“  -  2.15  x  10s  exp  M  0500/To ;  V-  2£/min) 


—  C*lculat#d 

•  Experimental 


io*s-k; 
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FIGURE  2  EFFECT  OF  SURFACE  REACTIVITY  ON  TEMPERATURE  PROFILE  IN  CATALYTIC 
MONOLITH  DUCT 

{Inlet:  gat  temperature  •  500  K,  0.5  vol.  %  C3Hgin  air;  k'  «  2.15  x  106  exp  (-10,500  TQ); 
k'  -  8.6  x  107  exp  (-8600/Ts):  V  -  28/minl 
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FIGURE  3  TEMPERATURE  PROFILE  IN  CATALYTIC  MONOLITH  DUCT 

Inlet:  G»s  temperature  ■  500K;  0.5  vol  %  CjHgin  air;  k'  *  8.87  x  103  exp 
(-8800/Ts);  V  ■  6£/min] 


INLET  SURFACE  TEMP.  PARAMETER/1.4  (Tso-500),  [K] 
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LOCATION  OF  TEMPERATURE  MAXIMUM  Xmax  IN  CATALYTIC  MONOLITH  DUCT 
AT  DIFFERENT  GAS  VELOCIT'SS 

[Inlet:  gas  temperature  ■  500  K;  0.5  vol.  %  C3Hg  in  air;  k“  ■  2. '.5  x  10s  exp  (-10,500/TQ); 
k*  -  8.67  x  103  exp  (-8500/TJ) 
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FIGURE  5  MAXIMUM  TEMPERATURE  IN  CATALYTIC  MONOLITH  DUCT  AT  VARIOUS  INLET 
SURFACE  TEMPERATURES  ANO  GAS  VELOCITIES 

(Inlat  gas  tamparatura  •  500  K;  0.5  vol.  %  C3Hgin  air;  k"  ■  2.15  x  106  axp  {-1Q,500/Tq); 
k”  -  8.87  x  103  axp  (-8600/Ts)] 


